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A stochastic background of gravitational waves from hybrid preheating 
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The process of reheating the universe after hybrid inflation is extremely violent. It proceeds 
through the nucleation and subsequent collision of large concentrations of energy density in bubble- 
like structures, which generate a significant fraction of energy in the form of gravitational waves. We 
study the power spectrum of the stochastic background of gravitational waves produced at reheating 
after hybrid inflation. We find that the amplitude could be significant for high-scale models, al- 
though the typical frequencies are well beyond what could be reached by planned gravitational wave 
observatories like LIGO, LISA or BBO. On the other hand, low-scale models could still produce a 
detectable stochastic background at frequencies accesible to those detectors. The discovery of such 
a background would open a new window into the very early universe. 



According to general relativity, the present universe 
should be permeated by a diffuse gravitational wave 
background (GWB) with a variety of origins, from un- 
resolved point sources (gravitational collapse of super- 
novae, neutron star and black hole coalescense, etc.) to 
relic stochastic backgrounds from early universe phase 
transitions, inflation, turbulent plasmas, cosmic strings, 
etc. [l|. These backgrounds have very different spectral 
shapes and amplitudes that may, in the future, allow 
gravitational wave observatories like LIGO, LISA, BBO 
or DECIGO [l[ to disentangle their origin. 

There are already a series of constraints on some of 
these backgrounds, the most stringent one coming from 
the large-scale polarization anisotropies in the Cosmic 
Microwave Background (CMB), which may soon be mea- 
sured by Planck, if the scale of inflation is sufficiently 
high [2|. There are also constraints coming from Big 
Bang nucleosynthesis [|[ , as well as from millisecond pul- 
sar timing [4j. Furthermore, it has recently been pro- 
posed a new constraint on primordial GW coming from 
CMB anisotropies ||. Most of these constraints come 
at very low frequencies (typically from 10 -18 Hz to 10~ 8 
Hz), while present GW detectors work at frequencies of 
order 1-100 Hz, and planned observatories will range from 
10" 3 Hz of LISA to 10 3 Hz of Advanced-LIGO [H], which 
could detect GW associated with early universe phenom- 
ena like first-order phase transitions [y, 0] , or cosmic tur- 
bulence [|[, if these occur around the electroweak scale. 

In this Letter we want to describe a new stochastic 
background of gravitational waves that may help open 
a window into the very early universe phenomena. Re- 
cent observations of the CMB anisotropies seem to sug- 
gest that something like inflation must have occurred in 
the early universe. The process by which the energy 
density driving inflation is converted into all the radi- 
ation and matter we observe today is called reheating, 
and corresponds to the true Big Bang of the Standard 
Cosmological Model. The first stage of conversion, pre- 
heating [ll|, is known to be explosive, and generates 
in less than a Hubble time the huge entropy measured 
today. In chaotic inflation, the coherent oscillations of 
the inflaton during preheating generates, via parametric 
resonance, a population of highly occupied modes that 



behave like waves of matter, which collide among them- 
selves and whose scattering leads to homogeneization and 
local thermal equilibrium. These collisions occur in a 
highly relativistic and very asymmetric way, being re- 
sponsible for the generation of a stochastic background 
of gravitational waves [13 , fl3| with a typical frequency 
today of the order of 10 7 — 10 9 Hz, corresponding to the 
present size of the causal horizon at the end of high-scale 
inflation. There is at present no chance to detect such a 
background, not even by resonant superconducting mi- 
crowave cavities [H[ • 

However, there are models like hybrid inflation in 
which the end of inflation is sudden [l5[ and the con- 
version into radiation occurs almost instantaneously. In- 
deed, since the work of Ref. [ll| we know that hybrid 
models preheat in an even more violent way than chaotic 
inflation models, via the spinodal instability of the sym- 
metry breaking field that triggers the end of inflation, 
irrespective of the couplings that this field may have to 
the rest of matter. Such a process is known as tachyonic 
preheating 0, [TlJ and could be responsible for copious 
production of dark matter particles [l8| , lepto and baryo- 
genesi s 1 1911 , topological defects [16| , primordial magnetic 
fields [2fj|. etc. Moreover, it was speculated in Ref. [2l| 
that in (low-scale) models of hybrid inflation it might be 
possible to generate a stochastic GWB in the frequency 
range accessible to present detectors, if the scale of infla- 
tion is as low as Hi n { ~ 1 TeV. However, the amplitude 
was estimated using the parametric resonance formalism 
of chaotic preheating, which may not be applicable in 
this case. In Ref. [l7| it was shown that the process of 
symmetry breaking proceeds via the nucleation of dense 
bubble- like structures moving at the speed of light, which 
collide and break up into smaller structures (see Figs. 7 
and 8 of Ref. 17]). We conjectured at that time that 
such collisions would be a very strong source of gravita- 
tional waves, analogous to the gravity wave production 
associated with strongly first order phase transitions [y] . 
As we will show in this Letter, this is indeed the case 
during preheating in hybrid inflation. 

Hybrid inflation models [l5| arise in theories of par- 
ticle physics with symmetry breaking fields ('Higgses') 
coupled to flat directions, and are present in many exten- 
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sions of the Standard Model, both in string theory and 
in supersymmetric theories Q. Inflation occurs along 
the lifted flat direction, satisfying the slow-roll condi- 
tions thanks to a large vacuum energy po- Inflation 
ends when the inflaton \ falls below a critical value and 
the symmetry breaking field cf> acquires a negative mass 
squared, which triggers the breaking of the symmetry 
and ends in the true vacuum, <j> = v, within a Hubble 
time. These models do not require small couplings in 
order to generate the observed CMB anisotropies; e.g. 
a working model with GUT scale symmetry breaking, 
v = 10~ 3 Mp, with a Higgs self-coupling A and a Higgs- 
inflaton coupling q gi ven by g = \/2X = 0.05, satisfies all 
CMB constraints [T(| , and predicts a tiny tensor contri- 
bution to the CMB polarization. The main advantage of 
hybrid models is that, while most chaotic inflation mod- 
els can only occur at high scales, with Planck scale values 
for the inflaton, and V^ 4 ~ 10 16 GeV, one can choose 
the scale of inflation in hybrid models to range from GUT 
scales all the way down to TeV scales, while the Higgs 
v.e.v. can range from Planck scale, v = Mp, to the Elec- 
troweak scale, v = 246 GeV, see Ref. HEG3]. 

Reheating in hybrid inflation goes through four well de- 
fined regimes: first, the exponential growth of long wave 
modes of the Higgs field via spinodal instability, which 
drives the explosive growth of all particles coupled to it, 
from scalars [l6| to gauge fields [l9[ and fermions [lH; 
second, the nucleation and collision of high density con- 
trast and highly relativistic bubble-like structures asso- 
ciated with the peaks of a Gaussian random field like 
the Higgs [l7|]; third, the turbulent regime that ensues 
after all these 'bubbles' have collided and the energy 
density in all fields cascades towards high momentum 
modes; finally, thermalization of all modes when local 
thermal and chemical equilibrium induces equipartition. 
The first three stages can be studied in detailed lattice 
simulations thanks to the semi-classical character of the 
process of preheating [22j , while the last stage is intrinsi- 
cally quantum and has never been studied in the lattice. 

In this Letter we use lattice simulations to study the 
generation of gravitational waves during preheating in 
hybrid inflation and analyse the dependence of the shape 
and amplitude of the spectrum of gravity waves on the 
scale of hybrid inflation, and more specifically on the 
v.e.v. of the Higgs triggering the end of inflation. Gravi- 
tational waves are represented by a tensor metric pertur- 
bation, g^ v = rj^ u + hfiu, in the transverse traceless (or 
radiation) gauge. Its equation of motion in this gauge 
is Oh^ = 16ttGT^, with the harmonic gauge condi- 
tion d^h^ — ensured by conservation of the energy- 
momentum tensor. In the radiation gauge we can fix 
^oo = 0, and the resulting field is the usual tensor 
gauge- invariant metric perturbation hij, which satisfies 
the evolution equation h'(j — V 2 /ijj = l6TrGTlij, with 
Hij the anisotropic stress tensor, sourced by both the 
inflaton and Higgs fields, Uy = Vj<^° Vj(f> a + Vi\ VjX — 
l/3^[(V</) a ) 2 + (Vx) 2 ]- We solve the evolution equa- 
tions of the gravity waves hij together with those of the 
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FIG. 1: The time evolution of the different types of energy (ki- 
netic, gradient, potential, anisotropic components and grav- 
itational waves for different lattices), normalized to the ini- 
tial vacuum energy, after hybrid inflation, for a model with 
v = 10 -3 Mp. One can clearly distinguish here three stages: 
tachyonic growth, bubble collisions and turbulence. 

other coupled scalar fields in a discretized lattice, assum- 
ing initial quantum fluctuations for all fields and only 
a zero mode for the inflaton, following the prescription 
adopted in Ref. [17]. We also included the GW backreac- 
tion on the scalar fields' evolution via the gradient terms 
/i y V i<f)V ' j4>, although for all practical purposes these are 
negligible throughout GW production. We then evaluate 
the mean field values, as well as the different energy com- 
ponents, see Fig.[T] For the energy in gravitational waves 

we use the expression (327rG) = (d^hj^ d^hl^r^j = 

| (d^hij d v h 1 ^ , where the expectation value is over a 
region sufficiently large to encompass enough physical 
curvature to have a gauge invariant measure of the GW 
energy [23} | , and we have expressed the average over the 
transverse traceless tensor hJ T in terms of the average 
over hij, the solution of the (traceless) tensor evolution 
equation. The fractional energy density in gravitational 
waves is then p gw /po = 4too/v 2 m 2 , which can be used 
to compute the corresponding density parameter today 
(with n Tad h 2 ~ 3.5 x 10" 5 ), 

f\w h 2 = f! rad h 2 1 (d hj? dohf^) , 

where we have assumed that all the vacuum energy p 
gets converted into radiation, an approximation which is 
always valid in generic hybrid inflation models with v <C 
Mp, and thus ff<m = yA v. We have shown in Fig. Q] 
the evolution in time of the fraction of energy density in 
GW. The first (tachyonic) stage is clearly visible, with a 
slope twice that of the anisotropic tensor Hy . Then there 
is a small plateau corresponding to the production of 
GW from bubble collisions; and finally there is the linear 
growth due to turbulence. Note that in the case that 
H <C m, the maximal production of GW occurs in less 
than a Hubble time, soon after symmetry breaking, while 
turbulence lasts several decades in time units of m . 
Therefore, we can safely ignore the dilution due to the 
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FIG. 2: We show here the comparison between the power 
spectrum of gravitational waves obtained with increasing lat- 
tice resolution, to prove the robustness of our method. The 
different realizations are characterized by the number of lat- 
tice points (N), the minimum lattice momentum (p m in) and 
the lattice spacing (ma). The growth is shown in steps of 
mAt = 1 for the lower spectra and mAt = 5 for the rest. 

Hubble expansion, until the universe finally reheats and 
the energy in gravitational waves redshifts like radiation 
thereafter. 

We then compute the power spectrum per logarithmic 
interval in GW by performing a Fourier transform of the 
energy density, fl gw — J df / f f2 gw (/), as a function of 
the frequency /, where f2 gw (fc) = k 3 p gw (k)/2ir 2 p c , with 
p c the critical density today. Since gravitational waves 
below Planck scale remain decoupled from the plasma 
immediately after production, we can evaluate the power 
spectrum today from that obtained at preheating by sim- 
ply converting the wavenumber k into frequency [121 ] , 
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We have shown in Fig. [5] the power spectrum of gravita- 
tional waves as a function of wavenumber k/m. We have 
used different lattices in order to have lattice artifacts 
under control, specially at late times and high wavenum- 
bers. We have checked that the power spectrum of GW 
follows (turbulent) scaling after mt ~ 40, and we can 
thus estimate the subsequent growth in energy density 
beyond our simulations. We have left for a future pub- 
lication [24[ the detailed analysis of turbulence in this 
system. 

We will now compare our numerical results with an- 
alytical estimates. The tachyonic growth is dominated 
by the faster than exponential growth of the Higgs 
modes towards the true vacuum [17|. The (traceless) 
anisotropic strees tensor grows rapidly to a value of 
order fc 2 |0| 2 ~ 10 -3 m 2 v 2 , which gives a tensor perturba- 
tion \hJ?h% T \ 1/2 ~ 16vrGw 2 (mAi) 2 10" 3 and an energy 
density in GW, p gw /p ~ M-kGv 2 (mAf) 2 lfr 6 - Gv 2 , 
for mAt ~ 16. In the case at hand, with v = 1CP 3 Mp, 
we find p gw /po ~ 10~ 6 at symmetry breaking, which 
coincides with the numerical simulations at that time, 
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FIG. 3: The sensitivity of the different gravitational wave 
experiments, present and future, compared with the possible 
stochastic backgrounds; we include the White Dwarf Binaries 
(WDB) [H and chaotic preheating (\<j> 4 ) [TJ] for comparison. 
Note the two well differentiated backgrounds from high-scale 
and low-scale hybrid inflation. 



see Fig. [T] The production of gravitational waves in the 
next stage proceeds through bubble collisions. Assum- 
ing the bubble walls contain most of the energy density, 
and since they travel close to the speed of light pTj . 
it is expected that the asymmetric collisions will copi- 
ously produce GW, like those of a strongly first order 
phase transition. In that quick estimate sug- 

gests that the fraction of energy density is given by [6j 
Pgw/po ~ 1/20(RH) 2 ~ 8tt/60 {RmfGv 2 - 2Gv 2 , of the 
same order or slightly larger than the previous stage, for 
the typical size of bubbles, R ~ 3to _1 , upon collision pTj . 
which again corresponds to what is observed in the nu- 
merical simulations, see Fig.[TJ The subsequent turbulent 
stage [13, [25[ is expected to further produce GW with a 
spectrum_that scales with time in a well defined manner, 
see also 
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fcl^Cfc) = q 2 Gv 2 T 1.0 k 2 exp( _ 25 k 2 T 

2ir z po 



where r = mt and p = 1/7 is the corresponding turbu- 
lent exponent [20L [25| . This spectrum has a maximum 
at k/m ~ 1, and falls as k 2 for small k until it reaches 
the maximum wavelength, k ~ H, corresponding to the 
minimum frequency today, / m i n ~ 5 x 10 10 Hz A 1 / 4 v /Mp. 
For the case we were considering in our numerical simu- 
lations, with v = 10 _3 Mp and A ~ g 2 ~ 0.1, we find the 
power spectrum of Fig. O 

We have plotted in Fig. [3] the sensitivity of planned 
GW interferometers like LIGO, LISA and BBO, together 
with the present bounds from CMB anisotropies (GUT 
inflation), from Big Bang Nucleosynthesis (BBN) and 
from milisecond pulsars (ms pulsar). Also shown are 
the expected stochastic backgrounds of chaotic inflation 
models like A0 4 [l2l [l3T | , as well as the predicted back- 
ground from two different hybrid inflation models, a high- 
scale model, with v = 10 _2 Mp and A ~ g 2 ~ 0.05, and a 
low-scale model, with v = 10~ 5 Mp and A ~ g 2 
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corresponding to a rate of expansion H ~ 100 GeV. The 
high-scale hybrid model produces typically as much grav- 
itational waves form preheating as the chaotic inflation 
models. The advantage of low-scale hybrid models of in- 
flation is that the background produced is within reach 
of future GW detectors like BBO 

To summarize, we have shown that hybrid models are 
very efficient generators of gravity waves at preheating, in 
three well defined stages, first via the tachyonic growth of 
Higgs modes, which act as sources of gravity waves; then 
via the collisions of highly relativistic bubble-like struc- 
tures with large amounts of energy density, and finally 
via the turbulent regime that drives the system towards 
thermalization. These waves remain decoupled since the 
moment of their production, and thus the predicted am- 
plitude and shape of the gravitational wave spectrum to- 
day can be used as a probe of the reheating period in the 
very early universe. The characteristic spectrum can be 
used to distinguish between this stochastic background 
and others, like those arising from NS-NS and BH-BH co- 
alescence, which are decreasing with frequency, or those 
arising from inflation, that are flat (27j |. 



For a high-scale model of inflation, we may never see 
the predicted GW background coming from preheating, 
in spite of its large amplitude, because it appears at very 
high frequencies, much beyond present experiments' sen- 
sitivities, where no detector has yet shown to be sensitive. 
On the other hand, if inflation occured at low scales, even 
though we will never have a chance to detect the GW pro- 
duced during inflation in the polarization anisotropies of 
the CMB, we do expect gravitational waves from pre- 
heating to contribute with an important background in 
sensitive detectors like BBO. The detection and charac- 
terization of such a GW background, coming from the 
complicated and mostly unknown epoch of rehating of 
the universe, may open a new window into the very early 
universe, while providing a new test on inflation. 
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